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Abstract 

We use a fluctuating hydro-dynamic approach to calculate the orientation fluctuations correlation 
functions of a thermotropic nematic liquid crystal in a nonequilibrium state induced by a stationay 
heat flux. Since in this nonequilibrium stationary state the hydrodynamic fluctuations evolve 
on three widely separated times scales, we use a time-scale perturbation procedure in order to 
partially diagonalize the hydrodynamic matrix. The wave number and frequency dependence of 
these orientation correlation functions is evaluated and their explicit functional form on position 
is also calculated analytically in and out of equilibrium. We show that for both states these 
correlactions are long-ranged. This result shows that indeed, even in equilibrium there is long- 
range orientational order in the nematic, consistently with the well known properties of these 
systems. We also calculate the dynamic structure of the fluid in both states for a geometry consistent 
with light scattering experiments experiments. We find that as with isotropic simple fluids, the 
external temperature gradient introduces an asymmetry in the spectrum shifting its maximum by 
an amount proportional to the magnitude of the gradient. This effect may be of the order of 7 
per cent. Also, the width at half height may decrease by a factor of about 10 per cent. Since to 
our knowledge there are no experimental results available in the literature to compare with, the 
predictions of our model calculation remains to be assessed. 
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I. INTRODUCTION 



Thermal fluctuations in an equilibrium isotropic fluid always give rise to short-range 
equal-time correlation functions, except close to a critical point. However, when external 
gradients are applied, equal-time correlation functions may develop long-range contribu- 
tions, whose nature is very different from those in equilibrium. For a variety of systems 
in nonequilibrium states it has been shown theoretically that the existence of the so called 

Hneric scale invariance is the origin of the long range nature of the correlation functions, 
, . For instance, for a simple fluid under a thermal gradient its structure factor, which 
determines the intensity of the Rayleigh scattering, diverges as q~ A for small values of the 
wave number q. This dependence amounts to an algebraic decay of the density-density 
correlation function, a feature that has been verified experimentally [3], j^J. 

Although studies of the behavior of fluctuations in nematic liquid crystals about nonequi- 
librium states are more scarce, some specific examples have been studied theoretically. Such 
is the case of the nonequilibrium situations generated by a static temperature gradient 5], 
a stationary shear flow j(| or by an externally imposed constant pressure gradient 8|, 
[jj. In the first two cases it was found that the nonequilibrium contributions to the corre- 
sponding light scattering spectrum were small, but in the case of a Poiseuille flow induced 
by an external pressure gradient the effect may be quite large. To our knowledge, however, 
at present there is no experimental confirmation of these effects, in spite of the fact that 
for nematics the scattered intensity is several orders of magnitude larger than for ordinary 
simple fluids lo| . 

In this work we present a model calculation based on a fluctuating hydrodynamic de- 
scription with a time-scale perturbation formalism to calculate analytically the equal-time 
correlation functions of the transverse and longitudinal orientation components of a ther- 
motropic nematic liquid crystal. We derive the explicitly space dependence of these quan- 
tities and show that they exhibit long-range order not only in equilibrium, but also in the 
nonequilibrium state induced when the liquid is subjected to a heat flux induced by a station- 
ary thermal gradient. We also evaluate the effect of the thermal gradient on the dynamic 
structure factor of the fluid. However, since to our knowledge there are no experimental 
results available in the literature for light scattering from a nematic in this steady states, 
to estimate the predictions of our model we used experimental parameter values similar to 
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FIG. 1: Schematic representation of a plane homeotropic cell with a constant thermal gradient 
along z direction. The inset shows the scattering geometry. The scattering angle is 9. 

those used in light scattering experiments for an isotropic simple fluid. In this way the effect 
of this long-range behavior on the dynamic structure factor of the fluid is examined for a 
geometry used in light scattering experiments. We find that the external gradient introduces 
an asymmetry of the spectrum shifting its maximum towards negative frequency intervals 
by an amount proportional to the magnitude of the gradient which may be of the order of 
~ 7%. Also, the width at half height may decrease by a factor of ~ 10%. 

II. MODEL AND BASIC EQUATIONS 

Consider a thermotropic nematic liquid crystal layer of thickness d confined between 
two parallel plates maintained at the uniform temperatures T\ and T2, respectively, in a 
homeotropic arrangement as depicted in Fig. 1. The transverse dimensions of the cell along 
the x and y directions are large compared to d. The hydrodynamic state of the nematic is 
specified by the velocity field, v (r ,t), the unit vector defining the local symmetry axis 
(director field), n(r ,t), the pressure, p(r ,t) and the temperature, T(r ,t). We assume 
that the nematic has reached a stationary state characterized by T ss , p ss and n ss = e z , 
without convection, e.g., T^ ss = 0. By symmetry p ss and T ss may only depend on z and 
if we neglect the variation of thermal conductivity with z, the stationary state is defined 
by the solution of the equations for p ss , p ss and T ss , which are obtained from the general 
nematodynamic equations for this geometry, see Appendix in Ref. [9]. As a result, a 
linear stationary temperature profile is established in the cell, T ss (z) = To (l + Bz^j, where 
T = (Ti + T 2 ) 1 2 and B = (3/d with (3 = d (dT ss /dz) /T . 
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We want to describe the dynamics of the spontaneous thermal (hydrodynamic) fluctua- 
tions, Sp (r , t) = p (r , t)—p ss , etc, around this stationary state. A complete set of stochastic 
equations for the space-time evolution of t 



general hydrodynamic equations 



nL 112J 



Q, 3 



;hese fluctuations is obtained by linearizing the 
and by using the fluctuating hydrodynamics 
formalism of Landau and Lifshitz |13j . In this work we will only calculate correlation func- 
tions in the bulk of the nematic, that is, for points ~r and r ' with z and z' far away from 
the boundaries, which means, \z — z'\ <C d and \z + z'\ <C d. Moreover, since the cartesian 
components of the director and velocity fields turn out to be strongly coupled, it will be con- 
venient to work in a representation where this coupling is minimum. This is accomplished 
by defining the transverse component of the director fluctuating field as its projection along 
the perpendicular direction to the k — h ss plane, i.e., 8h\ = k^h ss ■ ^ k x 5~n^j, where 
k± = [k 2 x + ky) 1 ^ 2 and the upper tilde denotes the space- Fourier transform of fields. Simi- 
larly, the longitudinal component, 5h 3 = k ■ Sit, is its projection in the k — h ss plane. 
In previous work we have shown that for a thermotropic nematic the modes associated 



with the director relaxation are much slower than the visco-heat and sound modes For 
typical material parameter values of a thermotropic nematic the ratio of these relaxation 
times is T orientation / T visco _ heat ~ 10 5 . Therefore, transverse and longitudinal director fluctua- 
tions relax to equilibrium much more slowly than fluctuations in temperature, pressure and 
velocity fields. Actually, the existence of widely separated time-scales may be exploited to 
eliminate the fast variables from the general dynamic equations obtaining a reduced descrip- 
tion in which only the slow variables are involved. The time scaling perturbation method 
introduced by Geigenmiiller et al. Q| may be implemented in order to diminish the 

couplings between nematodynamic fluctuations. It allows us to find, on the slow time-scales, 
a contracted description in terms of the slow variables only with a reduced dynamic matrix 
which can be constructed by a perturbation procedure. 

By using this method it can be shown that in the slow time-scale, director fluctuations 



5hi and 5h 3 obey the stochastic equations |X7| . 

dtSh^ = —oj^ ykj Sh^ — a^, \i = {1 transverse, 3 longitudinal (1) 
where we have introduced the following abbreviations 

uJi (t) = [K 2 k\ + K 3 k 2 z ) 



l 7 i(l + A) 2 fc 2 2 
+ 4 v 2 k\ + 



(2) 
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1 (Kk* I Kk>) fl I 1 7i[(l + AR 2 + (l-A)A:l] 2 



u 3 Ik 

The fluctuating sources are defined as 



k\ 



2 V 2 k\ + ^ 



(3) 



(4) 



478 = fc 



fcT +A;T | 1 (1 + A)^ + (1- A) fcj / 2 ~ _ ~ 



(5) 



where A = 71/72 is a non dissipative coefficient associated with director relaxation, ji and 
72 being orientational viscosities; i^i, K 2 and are the splay, twist and bend elastic con- 
stants, respectively; and u±, u 2 and u 3 are three shear viscosity coefficients. The stochastic 
components of the quasi-current of director field and the stress tensor, Tj and Ey, respec- 
tively, are zero averaged (Tj) = (Ey) = and in equilibrium they satisfy the fluctuation 
dissipation relations |9( 

1 



(Ti (T, t) Tj (7* t 1 )) = 2k B T-5±5 (r - 5 (t - t') 



7i 



(6) 



(Ey (~r, t) E im (~r', t')) = 2k B Ti> ijlm 5 (~r - ~r') 5 (t - t') , 
where k B is Boltzmann's constant, T is the equilibrium temperature of the nematic, 5. 



(7) 



Sij—nin s is a projection operator and i/j,-j m is the viscous tensor as given in Ref. 



Relations 



(jnj and (j7J) describe stationary Gaussian Markov processes. 

Now, we will assume that in the non-equilibrium steady state the fluctuation dissipation 
relations for Tj and E^ can be obtained from (JBJ) and (J7|) by replacing the equilibrium 
temperature T by the local ^-dependent stationary temperature T ss (z) . This assumption 
leads to 



Tj ( k , u ) Tj ( ', 



1 



7i 



<9 



2(27r) 4 fc B To-^ l + z£— U k + k')6(u + u'), (8) 



Ey ( fe , a; ) E /m ( fc a;' 



2 (2vr) 4 fc B To^ iJm ( 1 + iB^- ) 5 [ k + k ' ) 5 (u + u') . 



dk 



(9) 
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III. ORIENTATION AL CORRELATION FUNCTIONS 



The space-time Fourier transform of Eq. ((TJ is <5n M (jc, uo 

-l 



-iuj + ( k 



X m ( k , k';u,uj' 



and 



the director auto-correlation functions 



5h^ ( k',uj']} can be constructed in terms of the 



corresponding fluctuation-dissipation relations for a± and 03, 



X, 



er M ( A; , ) ( J 



/'/' 



-zcj + a;,,. ( k 



-iuj 1 + u;„ ( k ' 



(10) 



where no summation over the repeated index fi is implied. In these equations / denotes 
the space-time Fourier transform of the field /. The fluctuation-disspiation relations obeyed 
by a 1 and 03 may be found from Eqs. (JHJ) and Q. By inserting the resulting expressions 
into Eq. (|10p we obtain explicit forms for l w . First we shall examine the spatial limiting 
behavior of these correlations which is given by 



X^(T,T';t,t' 



-^y J dtdt'dudu'ei^^'^'-^-^X^ (y,~k';u,u'y (11) 



when (tf — t) — ► 0. Since the calculation of Xu and X33 is formally the same, we will 
only describe the procedure for X n . Calculating (&i ^ k ,u>j a\ yk ',u'j| with the help 
of Eqs. (jlj), (JHJ) and 0, replacing the result into (fTUj) and integrating over u', k ' and lu, we 
find that 

X n (r>, r"; t, f ) = X« (r , y ; t, t') + r"; t, . (12) 



where 



(i) 
11 



and 



V (2) 
A ll 



ik B T B 
'(2tt) 3 7i 



k B T ss (z) 
(2vr) 3 

71-^3 



c?A;- 



l^kT+K^ 



2 It - t'\ bi (k 



(13) 



3 iifc-(l ? -"r ?/ )-|t-t'|a;i(fc) 



Here we have considered t' > t and defined 

2 - - K2 



h k 



1 + A 



71^2^1 



#3 



1 l Tl (l + A) 2 ^ 
4 u 2 k 2 ± + v 3 k% 



(14) 



(15) 



2 7 fofci + ^ 2 2 ) 2 «i + ^ 3 ^ 2 

From Eqs. (|T3*j) and ((Hj) we can find the spatial range order of transverse director 
fluctuations by considering the limiting behavior at small times It — 1'\ or large distances 



\~r —~r'\, that is, when £ = K \ t — t'\/v\~r — ~r'\ « 1, which corresponds to the limit of 
static correlation functions. In this limit we find 



k B T ss (z) 

hm An = ; -rjo 

An(K 2 K 3 ) 1/2 



1 



r i — r 



7 |2 , K2 
■J- 1 + ^3 



1/2 



1 - 



z — z 



dz J 2T ss (z) 



, (16) 



where ("r^ — ~r' 1 _f' — (x — x') 2 + (y — y') 2 . 



Similarly, for the longitudinal director fluctuations we obtain 



lim \X 33 \ 



k B T ss (z) 



1 



(K./K,) 112 



An {K 1 - K 3 ) \\r - r ' 
' dT„\ z — z' 



r' ± \ 2 + ^(z-zr 



1/2 



X 



(17) 



dz J 2T ss (z)_ 

Note that in equilibrium, B = 0, T ss (z) = T, a constant, and the transverse orientational 



correlation functions reduces to 



lim \X°l\ = ^ 

1 An(K 2 K 3 ) 1/2 



?' ± \ 2 + %(z - z'f 



1/2' 



(18) 



while the corresponding correlation function for the longitudinal director components is 



\im\Xm 



knT 



4tt {K x - K 3 ) 




(K./K.f 2 



1/2 



(19) 



Note that both and A~|J are long ranged as could have been anticipated due to the 
well known property of a nematic which spontaneously exhibits a macroscopic orientational 
order. They decay algebraically as \~r — ~r'\ 1 . In Fig. 2 we plot the normalized static 
correlation in equilibrium X^ = And (K 2 K 3 ) 1 ^ 2 X^f/k B T for x — x' = y — y' = z = 0, T = T 
and as a function of normalized distance z'/d. In the stationary state, B ^ 0, the behavior 
of both correlations is modified by the presence of a term proportional to the temperature 
gradient which behaves as \~r — T*'|°, that is, which does not decay, in the direction of the 
temperature gradient. This is also depicted in Fig. 2, where we plot the normalized static 

1 /2 

correlation X u = And(K 2 K 3 ) 1 Xu/k B T for the same conditions as in the equiluibrium 
case and a value of the normalized thermal gradient (3' = dTL^ 1 (dT ss /dz) = 0.5 > 0. Note 
that the because of the presence of dT ss /dz, Xu becomes asymmetric. When compared to its 
equilibrium value, correlations increase in the direction of lower temperatures and decrease 
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FIG. 2: ( ) Decay of X n as a function of z'/d. ( ) Decay of In as function of z'/d for 

(3' = 0.5. The elastic constants values are Ki = 4.4 x 10~ 7 dyn , K% = 8.9 x 10~ 7 dyn. 

in the opposite direction. Indeed, for a plate separation d = 10~ 2 cm, the difference between 
both curves becomes significant for z'/d ~ 10" 1 . This means that the wavevectors sensitive 
to this difference are of the order of q ~ 10 3 cm -1 . Since for light scattering the wavevector k 
and the scattering angle 9 are related by q = 2qi sin 0/2, where ~ 10 5 cm _1 is the incident 
wave number, this implies very low scattering angles 9 ~ 0.1°. A quantitative evaluation of 
this nonequilibrium effect on a measurable property will be discussed in the next section for 
the structure factor of the fluid. 

IV. LIGHT SCATTERING SPECTRUM 

As an application of the previous theory we now calculate the light scattering spectrum 
of the nematic in the stationary state and for the scattering geometry defined in Fig. 1. For 
a nematic dielectric tensor fluctuations come mainly from director fluctuations and for the 
present model the spectral intensity of the scattered light is proportional to the dynamic 
structure factor which is given by the real part of (5hi (~q, uj) 5hi (— q , —to)), C (~q, cu), 



where V s and t s are the scattering volume and scattering time, respectively, q = q i — q s 
is the scattering vector (q y = 0) and uj = u s — u,i is the frequency shift. e a = e\\ — e± denotes 
the dielectric constant anisotropy. Evaluating this expression with the help of Eqs. (jHJ), © 



S(q,u) 




C(g>) 



(20) 
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we obtain 



_> ^_ 2elk B T cos 2 e g(Y) f 1 (dT ss \2uq^) 



7i ^ 2 + ( 5 ) \ T \ / cu 2 + u; 2 (9 
In this result we have considered the spectrum produced by a scattering volume located at 
the center of the cell. The nonequilibrium contribution has been written up to the smallest 
power of the wave number q, that is, up to the leading term in the hydrodynamic limit 
q — > 0, uji ( q ) is given by Eq. (J2J) and the functions a (q) and (3 (q) contain angular 
information through 

a{7) = 1+ (l±^V^- 2 (22) 
V 2 J v 2 ql + v z q 2 z 



and 



7i 



K 3 a(t) + {K 2 q 2 x + K 3 q 2 z ) 



2 



(z/ 2 g 2 + u 3 q 2 z ) 2 
In equilibrium, (j2"Tj) reduces to 



(23) 



2e 2 a k B Tcos 2 9 a(~q) 

7i ^ 2 + ^1 (7) 



g-n^)= a , ., (24) 



which behaves as q~ A . This dependence is responsible for the well known long-range order 
spatial behavior of the orientational correlations exhibited spontaneously by a nematic in 
equilibrium. On the other hand, the nonequilibrium contribution is 

S ne q ^ ^ = 2e l k B sin 2 6 f dT^X 2uq z at (~q) (~q) ^ 



7i \dz J [cu 2 + u 2 ("?)]' 

which also shows a long range-order decaying as q~ 5 . Note that close to equilibrium the 
size of the shift is indeed proportional to dT ss /dz. Its odd dependence on uj introduces an 
asymmetry in the shape of the structure factor, shifting the maximum towards the region 
of negative values of ou, as shown in Fig. 3. 

For a fixed scattering vector q , we define the dimensionless dynamic structure factor, 
So (uJq), in terms of the normalized thermal gradient j3' = dT^ 1 (dT ss /dz) and the normalized 
frequency u = u/u>i (q), by 

So M = f ( y\ = -Ml - P^s] (26) 

where A = q z (5 (~q) /oj\ (~q) d. In Fig. 3 we compare So (<^o) an d Sq 9 (uq) for low scattering 
angles, 9 ~ 0.1°, k{ ~ 10 5 cm -1 , f3' = 0.5 and typical values of material parameters of a 
thermotropic nematic |19]. From this result it follows that the increase of the maximum is 
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(B 

FIG. 3: ( ) Normalized structure factor at equilibrium Sq 9 (ujq). ( ) Nonequilibrium structure 

factor 5*o (wo)- We use f3' = 0.5, d = 10 cm, 6 ~ 0.1°, hi ~ 10 5 cm~ l . The values of the elastic 
constants are the same as in Fig. 2, 1/2 = 0.41 poise, ^3 = 0.24 poise, 71 = 1.03 poise, A = 1.03. 

about 7% and that the decrease in the half width at half height is 10%. This shows that 
the nonequilibrium state may induce changes in the dynamic structure factor which might 
be detected experimentally. 



V. DISCUSSION 



We have shown theoretically that the orientation correlation functions for a nematic 
liquid crystal exhibit long-range order both, in equilibrium and for a steday state induced 
by an external thermal gradient. We also estimated the influence of the stationary heat 
flux on the light scattering spectrum of a thermotropic nematic. The analysis carried out 
in this work included only the orientation correlation functions and the model has been 
constructed so that it corresponds to an experimental arrangement appropiated to detect 



the so called mode 2 of the spectrum 



10j . However, it should be emphasized that this is a 



model calculation and since to our knowledge there are no experimental results to compare 
with, the correctness of our choice of experimental parameters such as d, 9 or /3', remains to 
be assessed. However, as found in other nonequilibrium states for liquid crystals jjj, their 
magnitude suggests that they might be experimentally detected. 

It is worth pointing out that in previous work we have analyzed the effects produced by 
a nonequilibrium state induced by a concentration gradient in Ref. Other correlation 
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functions such as velocity-velocity or temperature-temperature may also be calculated by 
using the same approach. These correlations turn out to be also long-ranged in the presence 



of the thermal gradient as reported in Ref. |17(, |2fjj. They are also used to compare the 



physical mechanisms responsible for the appearing of long-range order in stationary states 
of nematics. 
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